Quasi-particle properties and Cooper pairing in trapped Fermi gases 
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The possibility for the particles in a Fermi gas to emit and reabsorb density and spin fluctuations 
gives rise to an effective mass and to a lifetime of the quasi-particles, as well as to an effective 
pairing interaction which affect in an important way the BCS critical temperature. We calculate 
these effects for a spherically symmetric trapped Fermi gas of ~ fOOO particles. The calculation 
provides insight on the many-body physics of finite Fermi gases and is closely related to similar 
problems recently considered in the case of atomic nuclei and neutron stars. 
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The experimental advances in cooling and manipula- 
tion of atomic Fermi gases are not only important in 
their own right, but also because of the possibility of 
providing a deeper understanding of the physics of finite 
many-body Fermi systems. A crucial feature of atomic 
Fermi gases is the possibility of exploiting Feshbach res- 
onances. These resonances occur when, due to the ac- 
tion of an external magnetic field, a two-body bound 
state crosses threshold. The result is a change in the 
sign of the two-body scattering length a and allows to 
explore Fermi superfluidity in all regimes from the weak- 
coupling limit, for a<0 and kp\a\ <C 1 (with kp being the 
Fermi wavenumber) to the strong-coupling one, for a<0 
and kp\a\ ~ 1, as well as Bose-Einstein condensation of 
two-fermion dimers (for a>0) 0. This crossover from 
fermionic to bosonic superfluidity has been the object 
of extensive work in the recent years and the first ex- 
perimental data are being collected 0. A fundamental 
theoretical problem is the evaluation of the critical tem- 
perature and the T — equation of state in the whole 
crossover region. In this paper we concentrate on the cal- 
culation of the effective quasi-particle interaction, mass 
and lifetime in a trapped Fermi gas as the resonance 
is approached from the fermionic, a < 0, side. As one 
moves closer to resonance the strength of the effective in- 
teraction increases and the quasi-particle properties are 
strongly renormalized due to the possibility of emitting 
and reabsorbing density and spin fluctuations. This prob- 
lem is particularly important when studying fermionic 
superfluidity, since both the superfluid critical temper- 
ature, T c , and the amplitude of the T = pairing gap 
depend strongly on the actual quasi-particle properties. 
These features are normally neglected in the standard 
crossover theories. As we shall show, they have a strong 
effect on T c . In what follows we develop a framework to 
take into account the effective quasi-particle properties 



and which is appropriate to finite size systems in a non- 
uniform isotropic trap. The case of uniform atomic Fermi 
gases has been dealt width in Ref. Q . In what follows we 
show that the discrete level structure has a crucial influ- 
ence on the problem, and point out the specific effects on 
T c of the various quasi-particle properties. We relate our 
approach to other calculations carried out for different 
Fermi systems. In fact the role of polarization effects on 
the quasi-particle properites has been recently studied in 
a number of fields of physics as nuclear physics, neutron 
stars and 3 He (cf. e.g. [j,|5| and references therein). 

We consider a Fermi gas of N particles in two internal 
states described by the Hamiltonian 

H = E / d 3 r ti(r)(H - ^ a (r) 

<T=Td 

+ gf d^$(r)$(r)^(r)$ T (r) (1) 

where H = -h 2 V 2 /2m a + V ext {r) with V ext {r) = 
m a uiQr 2 /2, loq being the trap frequency, and m a the 
atomic mass. We are here interested in the BCS side 
of the resonance where g < 0. We have furthermore as- 
sumed the number of particles and the external confining 
potentials V ex t to be equal for both internal states so that 
we also have /if = /if = /i. 

All the properties of the system are described in terms 
of the normal and anomalous propagators, which satisfy 
the Dyson equations (cf. e.g. ref. Q) 
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= ^ T T K) + ^ T T K)[s TT (w„)^ TT (w„) 
- w n K)jr/ T K)], (2 ) 

^ T K) = ^ i ° i (-^„)[VF i t T K)^ TT K) 

+ Sut-^j^K)], (3) 

where S%(w„) stands for Sf-|-f(x, x', cj„) and Sf-j-f (— o; n ) 
for 5fff(x',x, — u> n ), and similarly for the other propaga- 
tors. The quantity huj n = (2n + l)n/[3 is a fermionic 
Matsubara frequency with /3 _1 = fc^T, T being the 
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temperature of the system and ks the Boltzmann con- 
stant. S CT(T and Wf, - a are the normal and anomalous 
self-energies respectively, and the products of propaga- 
tors and self-energies imply integration over the internal 
spatial coordinates. Completely analogous equations to 
(J2J) and ® are satisfied by and =^tj.. The first step 
to be taken in carrying out the calculation, concerns the 
choice of a specific approximation for the self-energies. 
In this work we account for the possibility for particles 
to emit and reabsorb density and spin fluctuations and 
therefore assume T,^{uj n ) = gnj.(x)5(x — x') + T,^(oj n ). 
The first term corresponds to the usual w-independent 
Hartree self-energy. Treating the phononic contribution 
S|l to the self-energy in the Random Phase Approxima- 
tion (RPA) one finds: Eft(w n ) = -(ff7^)E m ^°k + 
w m )[n p (w m ) + 3n (Tz (a; m )-2no(a; m )]. Here n is the den- 
sity correlation function in the Hartree approximation, 
while n p and H az are the density and spin correlation 
functions respectively in the RPA. The relation between 
these functions follows introducing the Hartree correla- 
tion function for a single species Xoi to find Ho = 2xo to- 
gether with n p = 2Xo/[l-9Xo] and = 2xo/[l+ffXo]- 
Due to the isotropy of the system the spin-spin correla- 
tion functions in the x, y and z direction are equal, a fact 
which is accounted for by multiplying M„ z by a factor of 
3. Finally we also used ^ = Sf^ = The Hartree 
density correlation function admits the Lehmann repre- 
sentation 



n°(x,xV m ) = / — 



dJ 7T (X,X',U/) 



with 



7r°(x,x',a/) = £ 
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x (i\p{*)\j)(j\p{*!)\i)5{hw' -Ei+Ej). (5) 

In Eq. © p(x) = ^](x)^f (x) + -0](x)^(x), while i 
and j are the many-body eigenstates in the Hartree ap- 
proximation (particle- hole excitations), Ei and Ej are 
the corresponding energies and uj m is a bosonic Mat- 
subara frequency. n p admits a similar representation 
with i and j replaced by the many-body eigenstates in 
the RPA approximation (collective modes). The same 
is true for H az where in addition p must be replaced 
by a z = ■0|(x)-0|(x) — V?j(x)?/>j(x). The anomalous 
self-energy, on the other hand, is given by W<ti(w n ) = 
-/3" 1 E m ^ T f( Wm )^'TlK + w m ), where V^f \w m ) = 
gS(x — x') + (g 2 /4)[n p (a> r „) — 311^ (o; m )] is the effective 
quasi-particle interaction including the exchange of den- 
sity and spin modes. 

In a general spherically symmetric system it is conve- 
nient to express the Dyson equations in the Hartree-Fock 
(spherical) basis of eigenstates 0„(x) = R^f (r)Yi m (Q) . 



and magnetic quantum numbers) and Yi m is a spheri- 
cal harmonic. The corresponding energy measured from 
the Fermi level is denoted by = £ n ; independent 
of m. We shall also omit the spin indices. Keeping 
the matrix elements of the effective interaction only be- 
tween states connected by the time reversal V^,(w m ) = 
{fp\V eS (<jj m )\v' 'i?') and solving the Dyson equations in 
the proximity of T c one obtains 

W v {ui n ) = -J2 \ K e 5(^m)^'K + (6) 

V 1 771 

with #„(w n ) = W v (u n )G v (w n )G v (-u n ) and G~ x {uj n ) = 
ihcj n -^-T.P h {w n ), 

Introducing the Lehmann representation for the n's 
yields the following representation for the effective inter- 
action: 
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where g vu i — (vv\g8(-x. — yi')\v'v') is the ma- 
trix element of the bare interaction and v vv < (u) = 
(g 2 /4){vD\[ir p (pt,x.',w) - 37r CTa (x,x',o;)]|ivV). Corre- 
spondingly, after summation over the internal bosonic 
Matsubara frequencies, one obtains for the self-energy 

1 + n B (fru:) — n F {e v i) n B (hw) + n F (e v ') 



iuj„ — h 1 e,,> — uj 



iuj n — h 1 e v i + lu 



(8) 



with a vv i{u}) = (g 2 /A){iyiy'\[n p (-K 1 x',uj) + 37T CTai (x, x',w) - 
27r°(x, x', uj)]\v'i>), and ns(i) = (e l3x — l) -1 and np(x) = 
(e^+l) -1 being the Bose and Fermi thermal distribution 
functions respectively 0. Eqs. ©, Q) and JSJ consti- 
tute a set of equations of the type found in the theory of 
strong-coupling superconductivity (see for instance [a, 13) 
and should be solved self-consistently together with the 
number equation 



N = 



y f dx iimy^(x,x, Wn ) e ^', 



(9) 



which fixes the chemical potential. At low and interme- 
diate densities one expects the single-pole approximation 
for the propagators to be accurate. Introducing this ap- 
proximation for the spectral functions of G and & (see 
for instance 0,13) one obtains 
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i\huj n \ 



i\huj n \ 



(10) 



where e„ = + ReY, ph (u>)\f Uu=ev is the renormalized 



single-particle energy, 7„ = ImS p/l (w) 



the level 



Here v stands for n,l,m (radial, angular momentum width, and Z v = (1 — h 1 dYPj 1 (u>) / 'dio\% u) =cv) 1 < 1 the 
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quasi-particle strength. The single-pole approximation 
is accurate as long as the resulting single-particle level 
width is small compared with the Hartree level spacing 
and Z v is not too small. At these densities one also ex- 
pects the effect of the self-energies TP h and W on the 
number equation to be small compared with the Hartree 
one. We thus approximated Sf CTCT in Eq. © with the 
Hartree expression \(j)v(x)\ 2 / (ihu n — £„). Calcula- 
tions which go beyond the one pole approximation [Toj 
and which include vertex corrections 11] have been car- 
ried out recently to determine the pairing gap of super- 
fluid atomic nuclei. These calculations have shown that 
one recovers the results of the one-pole approximation for 
weak and intermediate coupling. 

Several further approximation schemes have been used 
in the literature. Introducing Eq. Q and (|1(J|) into Eq. 
© and neglecting the phonon contribution to the effec- 
tive interaction one recovers the approximations of Ref. 
|12|. For a Fermi gas this is in general a very poor ap- 
proximation since, for instance, in a uniform system the 
phononic contribution to the interaction produces a con- 
stant reduction factor (4e) 1//3 ~ 2.2 in the critical tem- 
perature in the weak-coupling limit |13|. On the other 
hand, if one keeps the phonon induced interaction and 
neglects the quasi-particle width one recovers the result 
of Ref. 0, found in the context of neutron matter. As 
we shall see, it is the quasi-particle width which has the 
largest effects on the critical temperature in the inter- 
mediate and strong coupling regimes. For these reasons 
we include both effects, together with the quasi-particle 
strength and the renormalization of the single particle 
energies, to get 
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FIG. 1: 7,, as a function of the Hartree quasi-particle energy 
£„ for A = 0.2 (solid line), A = 0.4 (dashed line) and A = 0.6 
(dot-dashed line). The quantities are measured in units of 
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with g vv i = Z v g vu i Z v i aixd A„ = Z u W v (e v 
sionless function h(e v t ,7^ ,T) is given by 



12] 



h(e u >,j u >,T) = 2k B T 
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The dimen- FIG. 2: T c barc and T c (in units of Hu>o) as a function of the 
interaction strength A. The values of T c Me are given in the 
caption of Table I, while those of T c include all contributions 
in Eq. ITU . 



(fiw„+7„') 2 + <v 



(12) 



and reduces to the usual hyperbolic tangent of BCS 
theory when 7^ = 0. The effective interaction is 
here treated in the Bloch-Horowitz approximation. This 
means that we have neglected the effect of the level width 
on the induced interaction and taken the absolute value 
of the external quasi-particle energy in the denominator. 
A detailed comparison between the results of this approx- 
imation and of the full self-consistent Gorkov theory was 
carried out in Ref. There it was shown that the two 
approaches lead to very similar results in the weak and in- 
termediate coupling regimes. Consequently we kept this 
level of approximation in our calculations. Eq. (|llf) is an 
eigenvalue equation for A„ and T c is the highest temper- 
ature at which the eigenvalue 1 appears. Notice that the 
equation contains a sum over v' which needs special care. 



The contribution of the bare interaction, in fact, leads to 
an ultra-violet divergence, since it arises from a contact 
approximation to the true interatomic potential. Several 
approaches have been developed to eliminate this diver- 
gence, approaches which involve a renormalization of the 
coupling constant g or the introduction of a pseudopo- 
tential 14]. For a trapped gas, for which only the bare 
interaction is considered, the results of these calculations 
may be reproduced by introducing a cut-off in the sum 
at e v i = ep ■ Because the main scope of the present 
work is on the phonon induced effects, we have treated 
the bare interaction at this level of approximation. Due 
to the explicit dependence on e u > , the phonon induced 
interaction has a natural cut-off. Numerical calculations 
have shown that the terms with \e v >\ J> £f give negligible 
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A=0.2 


A=0.4 


A=0.6 


V 


n no 


U.ou 


n on 

u.yu 


g+ ReE ph 


1.01 


1.03 


1.03 


9+ Z 


0.98 


0.90 


0.93 


9+ 7 


0.97 


0.85 


0.76 


V ca + RoE p ' 1 + 7 + Z 


0.88 


0.66 


0.61 



TABLE I: The table gives the ratio T , c /T c bare associated with 
the indicated contribution. J ,bare is calculated including only 
g and we find the following results: for A = 0.2, fcsT,? arc = 
5.03 x 10" 2 ftu , for A = 0.4, fc s T barc = 0.34/kj and for A = 
0.6, fc s T barc = 0.75ftwo. 



contributions. 

Applying the formalism to a spherically symmetric infi- 
nite square well we found the results of Gorkov and Melik- 
Barkhudarov [l3| in the weak-coupling limit. In the inter- 
mediate coupling regime we did not find the strong reduc- 
tion predicted in Ref . Q ; this may be due to the different 
set of approximations used. For a harmonically confined 
gas the finite level spacing strongly affects the properties 
of the system, and we find the results reported in Ta- 
ble I. We have performed the calculations for a cloud of 
~ 1000 particles and for three different coupling strengths 
A. Here A is defined in analogy with a uniform system 
as \g\m a k F /2TT 2 h 2 with k F = (Zn 3 ) 1 / 3 ^)} 1 / 3 , and n(0) 
being the central density of the cloud. 

In the weak-coupling regime the effects of the induced 
interaction and self-energy are strongly suppressed by the 
discrete shell structure, and the reduction predicted by 



Gorkov 13] for a uniform system is almost absent. This 
is because the largest contributions to the sum in Eq. 
(|llfl comes from the condition e„ = e u > = lu = 0, which 
corresponds to particles at the Fermi surface and to the 
exchange of a phonon with zero energy. The discrete 
level structure causes the values of the spectral functions 
relative to the II functions (see Eq. 10) to be negligible 
near ui = 0, and consequently the effect of the phonon- 
induced interaction to be small. 

The real part Re£ ph causes an increase in the crit- 
ical temperature since it leads to an effective increase 
of the density of levels at the Fermi energy, and thus 
to a stronger effective coupling strength. On the other 
hand, both Z v and 7^ depress the value of T c . The for- 
mer causes a weaker effective interaction between quasi- 
particles, as can be deduced from Eq. ltTT|) recalling 
that Z v < 1 0. Because 7^ measures the energy 
range over which the quasi-particle state is spread due 
to the coupling to vibrations (lifetime), its presence ef- 
fectively inhibits pairing between quasi-particles |12| . In 
Fig. n we show the value of j v for three different cou- 
pling strengths. Its effect is important in the interme- 
diate coupling regime, as can be seen in Table I for the 
cases A = 0.4 and A = 0.6. In Fig. [3]are shown the values 
of T^ arc , obtained when only the (bare) direct interaction 
g is included, and of T c when all the terms in Eq. i|Tl"|) 
are considered. 

We acknowledge useful discussions with P.F. Bortignon 
and the technical support of F. Marini. 
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